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Abstract 

We study theoretically the transition temperature T C (H) of a qu asi- 
one-dimensional anisotropic superconductor in a magnetic field, which is 
thought to be realized in organic conductors, (TMTSF^X. In the weak 
field limit, — dH C 2/dT for the anisotropic singlet is shown to be V3 times 
larger than that for the anisotropic triplet in the case of attractive in- 
teraction of electrons at nearest-sites along the c axis and magnetic field 
along the b axis. In the strong field region the transition temperature 
for spin-triplet is restored to zero field critical temperature while the an 
isotropic spin-singlet is suppressed by the Zeeman effect. There should be 
tra nsition from anisotropic singlet to anisotropic triplet as the magnetic 
field is increased. 

1 Introduction 

Recently, the reentrance of the superconductivity of quasi-one-dimensional elec- 
trons in a strong magnetic field attracts theoretical interest, [[l] g, ||, [| The 
anomaly of the resistivity in a strong field is observed, || which is thought to 
be a signal of the superconductivity. The superconductivity in a strong field 
is understood as follows. When the magnetic field is applied along the b axis, 
which is different direction from the field-induced spin density wave (FISDW) 
case, the semi-classical orbits of electrons are localized in the a-b plane (where 
a and b is the most conducting and second conducting axes, respectively), i.e. 
the motion along the c axis is bounded. Since the orbital frustration due to 
the magnetic field comes from the motion of electrons in the plane perpendic- 
ular to the magnetic field, electrons can make Cooper pairs without affected 
by the orbital frustration in the strong field limit. The superconductivity in 
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strong magnetic field can be also understood by nothing that the quasi-one- 
dimcnsional system in the magnetic field is essentially one -dimension with the 
field-dependent effective interaction. [p| 

Lebed' Q and Dupuis et al. |^|. ||,]4|] have calculated the mean-field critical 
temperature of a quasi-one-dimensional superconductor. In these papers they 
considered both singlet and triplet pairings but they assumed the isotropic su- 
perconductivity, i.e. the energy gap is constant at the Fermi surface. However, 
considering the recent experimental results of NMR relaxation rate Q, impu- 
rity effect § and proximity to the SDW phase § in the (TMTSF) 2 X, the 
superconductivity in the quasi-one-dimensional organics is thought to be not 
the conventional type with the constant energy gap. 

In the previous paper we have calculated the effective interactions in the 
magnetic field by assuming the attractive interaction between electrons at the 
nearest-sites in the c direction, which results in the anisotropic superconductiv- 
ity. H For this case the transition temperatures for the spin-singlet and spin- 
triplet are degenerate in the absence of the magnetic field. We got the different 
magnetic field dependence of the superconductivity of anisotropic singlet and 
triplet. In that paper we took the quantum limit approximation (QLA), which 
is valid for the strong field. In the intermediate field, the transition temperature 
calculated in the QLA is known to be smaller than that without the QLA. [jlO| 

In this paper we study the transition temperature of the anisotropic su- 
perconductivity in a quasi-one-dimensional system by extending an elegant ap- 
proach used by Lebed' and Dupuis et al. 

We also consider the effect of Pauli pair-breaking for an anisotropic super- 
conductivity by taking account of the Zeeman term. In general, the critical 
field of spin-singlet superconductivity is decreased by the Pauli pair-breaking 
effect. When the number of quasi particle polarized by magnetic field exceeds 
a critical value, Cooper pairs (k + q/2 f, — k + q/2 J.) with total momentum 
q are stabilized rather than those of (k |, — k J,). The inhomogeneous super- 
conductivity with spatial oscillation of the gap function is then possible, as was 
pointed out by Fulde and Ferrell and Larkin and Ovchinnikov (FFLO 
state). In three dimensional case with a spherical Fermi surface, the FFLO state 
is known to be possible only in the small region of magnetic field. However, in 
the one-dimensional case, the electron (— k + q/2 |) is always on the down-spins 
Fermi surface for any (k + q/ 2 \) on the up-spin Fermi surface if we take an 
appropriate choice of q. (l3[ Q As a result, at least a half of the Fermi surface is 
available for pairing, so that the reentrance in the strong field is not completely 
suppressed by the Pauli pair-breaking effect even for the spin-singlet pairing. 
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2 Gap equation 



We consider the anisotropic superconductivity described by the dispersion law 
(we take U = k B = 1 throughout the paper) 

Eka = v F (\k x \ - ftp) - 2tbCOs(bk y ) - 2t c cos(ck z ) + a^ B H, (1) 

where vp = 2at a s'm(akF) is the Fermi velocity for the motion along the a axis, 
t a , tb and t c are the transfer integrals along a, b and c axes, respectively, a, 
b and c are lattice constants and a\x B H is the Zeeman energy for \ (J.) spin 
(a = +(—)). For simplicity we take a, b and c axes to be perpendicular each 
other and to be along the x, y and z directions, respectively. 

For the non-interaction Hamiltonian Ho of a quasi-one-dimensional system 
we take account of the magnetic field by the Peierls substitution, i.e. Hq = 
E(k — > — iV — eA). When the vector potential A is taken as A = (0, 0, —Hx), 
the eigenstates and the corresponding eigenvalues are given by 

{x I TO.) = e ik * x + ik y bl + ik ' cm - ia izf sin(k z c-Gx) ^) 

and 

e k, CT = v F{ak x - k F ) - 2t b cos(k y b) + a^i B H, (3) 

where the integers / and to label the cites in the directions y and z, a = sgn(fc x ) 
refers to the right/left sheet of the Fermi surface, and co c = evpcH. The one- 
particle Green's function of a quasi-one-dimcnsional system is represented by 

G%(x,x',ky,k z ,U> n ) = e -^[^(^o-G X )- S in(k z c-G X ')] 

x J2e ik ^-^G^(k x ,k y ,u; n ) (4) 



with 

G"(k x ,k y ,U n ) = - -. : : . j- rr — (5) 

vjj n — VF(ak x — kp) + 2% cos(k y b) — o\i B H 

where u n = (2n + l)nT is a Matsubara frequency. 
The interaction Hamiltonian is written as 

U ^ = \ E E ^CVi)€<Vi)e(r,)C(r;), (6) 

a,a' ,a.<7' <i,j> 

where Uij is the interaction between electrons at i and j sites and a — — a. The 
^>"(r)'s are operators for electrons moving on the sheet a of the Fermi surface: 

CM -EE ^k»C£ ;<T , (7) 

Q^a; >0 ky,k z 
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where C£ a (C^ a ) is the annihilation (creation) operator of an electron in the 
state 4>^Jr). 

In this paper we take the nearest-site interaction along c axis and on-site 
interaction as 

{U if Ti = Tj ± cz 
U ifr 4 =r, (8) 
otherwise 

We define the order parameter corresponding to spin-singlet pairing state 
A; s i ) (r i ,r j ) = [/ i ^(Vf(r j )Vr(^)) (9) 

a 

and spin-triplet pairing state 

A« ( Ti , rj ) = Ui^WZiTjWWa. (10) 

a 

In the mean field approximation, the linearized gap equation is obtained as 
A(x,q z )= f dx'K(x,x',q z )A(x',q z ), (11) 

J\x-x'\>d 

where <i _1 is related to cutoff energy 0. The kernel K for the anisotropic spin- 
singlet (A = Af jj is given by 

COS {CK Z )*~*^ \X , X , fcyi f^z 

w n a,k y ,k z 

x G| (x,x , ky,—k z H — , w n ) 

-f7T cos[2u B iiXx- x')/«f] 
_ 6cw| sinh[|a; - x'^ttT/wf] 

x [Jo (A(ar, x', q z )) + J 2 (A(x, x 1 , q z ))] , (12) 

where J n {A) is the nth-order Bessel function, G = eHc, and 

A(x,x',q z ) = — sin[^-(x- x')]sm[q z ^- - — (x + x')]. (13) 
For the anisotropic spin-triplet (A = or An) we get the kernel as 

K^(x,x',q z ) = -2UTY, E s™ 2 {ck z )G a a (x, x', k y , k z + |, u n ) 



K^(x,x',q z ) = -2UTJ2 E ™s 2 (ck z )G^ (x , x' , k y , k z + ^, u n ) 



<^n. a,ky,k z 

Qz 
2 

-C/T 1 

feci;! sinh[|a;-x'|27rT/uF] 

x [J (A(x,a;', gz )) - J 2 {A(x , x' , q z ))] , (14) 
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while for the isotropic spin-singlet the kernel is given by jl], ^, || 

K^Hx x' a ) - - U o Tcos i 2 ^H(x~x')/v F ] 

[X,X,qz) ~ bcv* S mh[\x-x>\2nT/v F ] Jo[A[X,X,qz)) - ( &j 

If the on-site interaction is repulsive (Uq > 0), the isotropic spin-singlet is not 
realized. 

In the present assumption for the interaction the isotropic spin-triplet does 
not occur. The isotropic spin-triplet is possible only when the backward and 
forward scattering constants are different, which will be the case when the at- 
tractive interaction exists between electrons on the nearest-sites along the a 
axis. 

2.1 Transition temperature for a weak field 

In the weak field limit (oj c <C T), where the magnetic length 2ir/G is much 
larger than the thermal length v f /2t:T, we can expand the order parameter 
and Bessel function in (x — x'). We can also neglect the Zeeman term in this 
limit. Then the second-order differential equation is given by 



v 



|^ + 4^[1 - cosfec - 2Gx))A = i^T c 2 (l - ^ ) A. (Hi) 



where T c q — (vF/ird) exp(— irbcvF/U) is the transition temperature in the ab- 
sence of magnetic field, Q{z) is the Riemann's zeta function and 

1 for isotropic superconductivity 

i for anisotropic spin-singlet (17) 
| for anisotropic spin-triplet. 

When T c q <C t c , it can be further simplified by replacing the periodic potential 
in the lhs by a set of decoupled harmonic potential located at points (q z c/2 + 
mr)/G. Then we get the critical temperature as 

rm^-^m^H, (is) 

or 

AH b c2 _ 1 8^ 2 T c0 

dT y/2!j7({3)t c VFec' [ 1 

The upper critical field H h c2 in this region is \[Z times larger for anisotropic 
spin-singlet than for anisotropic spin-triplet. If the magnetic field is ap- 
plied along the c axis, -AH c c2 /AT = 8tt 2 T c0 / '(7V2((3)t b v F eb), for both 
anisotropic spin-singlet and spin-triplet in the weak field limit. Its value 
is the same as that for the isotropic supeconductivity. In contrast with 
this case, if the attractive interaction works between electrons along b axis, 
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—dH£ 2 /dT — 87r 2 T c o/(7\/2rj(3)tbVFeb) for the weak magnetic field along c axis 
and -dH b c2 /dT = 8Tr 2 T c0 /(7V2({3)t c v F ec) for the field along b axis. 

We compare above results with the upper critical field (-diJ c b 2 /dT = 
1.84 and -dH^/dT = 0.111 [T/K]) observed by Murata et al p| for 
(TMTSF) 2 C10 4 . When we take the parameters as b = 7. 7 A, c = 13.5A, the 
ratio of transfer integrals are t c /tf, = 0.034 if we fit the experiments with the 
isotropic superconductivity. If we assume the attractive interaction along the 
c axis as studied in this paper we get t c /t b = 0.048. In the case of attractive 
interaction along the b axis, we get t c /tb = 0.028 to fit the experiment. 



2.2 Transition temperature for any magnetic field 

Without using the weak field approximation studied in the previous section, we 
can solve the linearized gap equation eq.(ll) numerically as done by Dupuis and 
Montambaux 0, [| for isotropic superconductivity. We write 



Ag(a:) 



A?,e i2,Gl , 



(20) 



where Bloch wave vector Q is taken as — G < Q < G. Then eq.(ll) is written 
as a matrix equation 



a<2 — aQ aQ 

^2f " / t ^ 1 2i.2i'^2i" 



where 



A 2l,2V = J2 K N -2l,N-2l>K(Q + NG), 
N 

In the above K{q x ) is given by 

K(q x ) = ^2 G^(k x ,k y ,uj n )Gg(q x - k x ,-k y ,-uj ri 



(21) 
(22) 



zN(0) 



E 



Re* I - 



In 



•Ml 



2^7 
~7rT~ 

4iirT 



(23) 



for the spin-singlet, where N(0) = 1/irvpbc is the density of states for one spin 
at the Fermi energy, 7 is the exponential of the Euler constant, f2 is the cut-off 
energy, and * is the digamma function. For the spin-triplet, Pauli pair breaking 
term 2(IbH in the digamma function should be eliminated. 

The coefficients Kn u n 2 are given for the anisotropic singlet as 



G 



<U - ^' il "'- , J (*£ C082( » )J -(^ 8i " ( » ) ) 

x J Nz ( — sin(u) ] , (24) 



and for the anisotropic triplet as 



x Jjv 2 ( — sin(u) 1 , (25) 



while for the isotropic singlet 



,(iso) -U Q hl {Ni - N2) f 27r du (tt c , \ (At c 



If we take the most divergent term (vf(Q + NG) + = 0) and neglect the 

others in eq.(22), we get the result in QLA. 

A. Without Zeeman effect 

We first consider the case without Pauli pair breaking. In Fig.l the transition 
temperature for anisotropic and isotropic spin-singlet is shown as a function of 
h = H/Hq, where Hq = 4t c /v-pebc. jl6| We take parameters for the organic 
conductor as a k 7.3l, c » 13.5A, t a » 2900K, t c w 17K and ak F = 7r/4, 
and get Ho « 11. 4T. Here two transition lines are characterized by the order 
parameter for Q — or Q — G. The very small oscillation of the transition line 
for the anisotropic spin-singlet exists in 0.05 < h < 0.25 in the present choice of 
parameters. The recovering to zero field critical temperature starts from lower 
field in the anisotropic spin-singlet case than in the isotropic case. In Fig. 2 
we plot the transition temperature for anisotropic and isotropic triplet with 
the same parameters as in Fig.l. The magnetic field dependence of transition 
temperature for isotropic spin-triplet is the same as that for isotropic spin- 
singlet. In the case for the anisotropic spin-triplet, the oscillation is large and 
the recovering starts from higher field than the isotropic type (Fig. 2). 
These results are consistent with those obtained in QLA. 



B. Zeeman effect 

We next consider the effect of Pauli pair-breaking by taking account of the Zee- 
man term. As seen in eq.(23), logarithmic divergences occur at Q = ±.2^11 /v-p 
and Q = ±(G — 2/ibH/vf) even in the presence of the Zeeman term. We take 
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these values of Q and calculate the transition temperature. As shown in Fig. 3, 
the transition temperature of the spin-singlet is considerably reduced by the 
Pauli pair-breaking effect. The transition temperature for Q = ±(G—2iibH/vf) 
is too small to be seen for h > 0.1 and almost degenerate with that for 
Q = ±2^B-ff/fF for h < 0.1 in Fig. 3. In the case of equal-spin-pairing A^j 
and A||, Zeeman splitting does not change the transition temperature. Thus, 
one will observe spin-triplet superconductivity in the strong field. 

Even with the Zeeman effect, Cooper pairs of whole electrons in the Fermi 
surface contribute to superconductivity if the N("f )th energy level coincides with 
the N+l(|)th level by Zeeman splitting. In this case the Pauli pair breaking 
effect becomes not important, as shown in Fig. 4. This is the case when 2/iB-ff = 
Gvp as seen in eq.(23), which yields act a — ^o/ib/v^tt, where </>o = hco/e is 
the flux quantum. This condition is satisfied if act a is about 1/10 times of 
its value for organic conductors. This may be possible for a system with a 
large effective mass, or in the tilted magnetic field in the a-b plane, because the 
Zeeman splitting is in proportion to the magnetic field, while magnetic length 
with respect to orbital motion is in proportion to perpendicular component of 
the magnetic field. On the other hand, in the case of tilted magnetic field in the 
b-c plane, since hopping along b axis is suppressed at the same time, one has to 
consider competition between superconductivity and FISDW. 

In the weak field limit the field dependence of the critical temperature is 
given by the parameter t c v-p oc t c t ai while in the strong field region it is given 
by Hq oc t c /vF oc t c /t a and T c q. In order to fit the observed initial slope 
— dH^ 2 /dT\T c0 ||, we take t a t c ~ 2100K 2 for the anisotropic singlet case and 
t a t c ~ 1500K 2 for the isotropic singlet case. In Fig. 5 we plot the transition tem- 
perature for the isotropic singlet and isotropic triplet by taking t c = 1.5K and 
T c o = 1.36K. If we assume that the transition temperature for isotropic triplet in 
the absence of magnetic field is much smaller than that for the isotropic singlet, 
the result seems to be consistent with experiment as discussed by Dupuis j| 
and Lee et al. J^] In Fig. 6 we plot the critical temperature for the anisotropic 
singlet and anisotropic triplet as a function of the magnetic field by taking 
t c = 1.5K and T c q = 1.36K. In this choice of parameters oscillation of the 
transition temperature for anisotropic triplet is pushed to the lower field. The 
transition temperature for the anisotropic spin-singlet with Q = G — 2/ibH/vp 
changes discontinuously to very small value at H = H c ~ 3. IT. By changing Q 
we may get the continuous transition temperature in H > H c , but this will be 
smaller than the transition temperature for Q — 2[1bH j 'v-p. Thus we did not 
consider the change of Q. As shown in Fig. 6, the transition from anisotropic 
spin-singlet to anisotropic spin-triplet occurs about H — 1.5T and the transi- 
tion temperature increases above this field. This feature is not observed in the 
experiment, [js) This may be due to the assumption of the attractive interaction 
we have used in this paper. The part of attractive interaction would be caused 
by the anti-ferromagnetic spin fluctuation, which results in the spin- dependent 
effective interaction on the nearest-sites and favors the anisotropic singlet pair- 
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ing rather than anisotropic triplet pairing. Then the transition temperature for 
the anisotropic triplet becomes lower and the transition from anisotropic sin- 
glet to anisotropic triplet will occur at higher magnetic field. Then we can also 
explain the experiment by Lee et al. with the anisotropic superconductivity. 

3 Conclusion 

In this paper we have studied the transition temperature for a quasi-one- dimen- 
sional anisotropic superconductivity in a magnetic field applied along the b axis. 
Here we have assumed the nearest-site attractive interaction between electrons 
along c axis, which is independent of spin. The transition temperature for the 
anisotropic spin-singlet is same as that for the anisotropic spin-triplet in the ab- 
sence of magnetic field. We have shown that the critical field for the anisotropic 
spin-singlet is V3 times larger than that for the anisotropic spin-triplet in the 
weak magnetic region. 

We have calculated the transition line of spin-singlet superconductivity con- 
sidering the Zeeman effect. Although spin-singlet sate exists even with the 
Zeeman effect, the transition temperature is suppressed. On the other hand, 
spin-triplet superconductivity is not affected by the Zeeman effect. Therefore, 
the transition from spin-singlet to spin-triplet will occur as the magnetic field 
is increased. 

The present results for anisotropic spin-singlet are consistent with the recent 
experiment by Lee et al j5j, which shows the signal for high field supercon- 
ductivity, if we take into account the attractive interaction for the anisotropic 
spin-singlet and the relative suppression of the anisotropic spin-triplet caused 
by the anti-ferromagnetic spin fluctuation. Even in this case the transition 
from anisotropic spin-singlet to anisotropic spin-triplet will occur in the strong 
magnetic field as long as the spin-triplet is not destroyed completely. 
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Figures 



Fig.l. Critical temperature T c (H)/T c q vs magnetic field along the b axis for the 
anisotropic singlet (solid lines) and isotropic singlet (broken lines). The Zeeman 
effect is neglected. 

Fig. 2. Critical temperature T C (H)/T c0 vs magnetic field for the anisotropic 
triplet (solid lines) and the isotropic singlet (broken lines). 
Fig. 3. The same as Fig.l with the Zeeman effect. 

Fig. 4. The solid lines show Energy spectrum of a quasi-one-dimensional 
system in a magnetic field and the broken lines show the Zeeman splitting. 

Fig. 5. Critical temperature vs magnetic field for the isotropic singlet (solid 
lines) and isotropic triplet (broken lines). Parameters are taken to fit the initial 
slope of the critical field —dH^ 2 /dT to the experiment and different from that 
in Figs. 1^3. 

Fig. 6. Critical temperature T C (H)/T c0 vs magnetic field for the anisotropic 
spin-singlet (solid lines) and for the anisotropic spin-triplet (broken lines) with 
the different choice of parameters from that in Figs. 1^3. The Zeeman effect is 
taken into account. 
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